Spitzer IRS observations of k+a galaxies: A link between PAH 
emission properties and AGN feedback? 
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ABSTRACT 

We have performed IRS low resolution 5-12 micron spectroscopy on a sample 
of galaxies selected to be at 3 distinct post-starburst evolutionary stages based 
on their optical spectral indices. The resulting IRS spectra show distinctive PAH 
emission line structures at 6.2, 7.7, 8.6 and 11.3//m and little sihcate absorp- 
tion, indicative of ongoing star formation. However the PAH inter-line ratios, 
in particular the 11.3/6.2/im and 7.7/6.2/im ratio, show large variations. These 
variations are found to correlate with both time since the most recent starburst 
and AGN activity. We speculate that the evolution observed in these PAH ratios 
is related to an increase in AGN activity with time since star burst. 

Subject headings: galaxies: evolution — galaxies: starburst — infrared: galaxies 



Introduction 



The nature of the connection between starbursts and AGN has been long debated. 
While the effectiveness of AGN feedback in quenching star formation has been heralded 



as a key component of ga laxy formation by semi-analytic models (e.g. ICroton et al.l 12006 



De Lucia &: BlaizotI 120071 ) . it is only very recently that observation al evidence of AGN ac- 
tivity co-incident with post-starburst signatures has b een found (ISchawinski et al.l 120071 : 
Bundy et al.ll2008l : lKaviraill2008l : ISchawinski et al.ll2008l ). However these discoveries form far 
from a consensus, with some suggestions that AG N feedback alone niay not be suffi cient for 
the truncation of star formation in these galaxies (IBundy et al.ll2008l : iKavirajl 120081 ). 



The mid-IR contains unique diagnostics of the star formation, AGN activity, and dust 
obscuration levels of galaxies. Rest-frame 5 — 12//m observations can decipher the relative 
strength of these attributes via investigation of the PAH emission properties, continuum 
slope, and silicate absorption strength, respectively. 
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Here we utilise Spitzer InfiaRed Spectrograph (IRS) observations of a sample of galaxies 
which appear at different stages of post-starburst evolution as determined by their optical 
spectra. In particular we use the B.5, [OII]3272A and 4000A break spectral features to iden- 
tify k+a galaxies at 3 different epochs, k+a galaxies are typically identified via their optical 
spectrum, with strong absorption in the Balmer series lines, but no significant emission in 
the [Oil] 3727 A line or other emission lines associated with star formation (Dressier & Gunn 
1983; Couch and Sharpies 1987; among others). This kind of optical spectrum is indicative 
of a galaxy which has undergone a very recent burst of star formation as the young, short- 
lived (< 1 — 2 Gyr) A and B stars are still observed, but no ongoing, or very low level, star 
formation. Spectral modelling of k+a galaxies typically suggests a star-burst in which 1-10% 
of the galactic stellar mass is in the form of new stars (i.e. Balogh et al 1999). 



Data 



We selected 12 galaxies in three categories from the overlap between the Spitzer SWIRE 
survey and the SDSS spectroscopic dataset for observations with IRS. The targets were se- 
lected via their [OII]3727A, R5 (410lA), and D4000 indices to be at distinct post-starburst 
stages. The B.S and D4000 indices are well established as good indicators of stellar age (i.e. 
Balogh et al]ll999l : iKauffmann et al1l2003l ). while the [OII]3727A is a well known indicator 
of current star formation. Measurements of the indices were performed as per Roseboom 
et al. (2006). To determine our s electio n criteria we turn to the stellar population synthe- 
sis models of iBruzual &: CharlotI (120031 ). Models which assume an exponentially declining 
starburst (r = 100 Myr) on top of a passively evo lving single stellar population of 11 Gyr 
were generated using the lBruzual fc CharlotI (l2003l ) code. Burst strengths of 1,5 & 10% were 
considered. From these models we determine 3 sets of & D4000 criteria which isolate 
galaxies with roughly 0.03, 0.3 and 1 Gyr of post-starburst evolution. In the case of the 
later two selections an additional [Oil] selection cut is imposed to filter out galaxies with 
significant ongoing star formation. The specifics of this selection are; 



Line Young k+a k+a Old k+a 

[OII]3727A(A)* - <6 <6 

H5 (A) > 4 > 5 > 2 

D4000 1 - 1.2 1.3 - 1.5 1.6 - 1.8 

Time since burst (Myr) 25 300 1000 
*In emission 
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Observations were performed with the short-low module in staring mode, two nod po- 
sitions per spectral order with an offset of 20" were observed with ramp times between 
400-600s. The data were reduced using a combination of the SMART (Hidgon et al. 2004) 
and SPICE software environments. First individual ramps of the same order and nod were 
combined in SMART utilising a median filtering algorithm. Sky subtraction is performed 
via subtraction of different order observations at the same nod position. Once combined 
and corrected for background/sky contributions the spectra are extracted using the SPICE 
pipeline assuming these objects are point sources for Spitzer at these wavelengths. Finally 
the two nods are combined in SMART to produce the final spectrum. 

11 of the 12 IRS spectra are found to have acceptable S/N spectra (> 5) however 
one of the IRS observations appears to have featureless, flat spectrum, suggesting that the 
exposure times were not sufficient to detect the mid-IR features, or the IRS slit was misplaced 
in relation to the mid-IR emission. This failed observations is omitted in the following 
discussion. 



Results 



Figure [H shows the IRS spectra for the 11 sources discussed here. For comparison the 
starburst template of Brandl et al. (2006) is also shown. All the spectra, with the exception 
of Oka 1, show clear PAH emission features at 6.2, 7.7 and 11.3/im. Oka 1 is found to have a 
completely stellar dominated thermal spectrum, similar to early-type galaxies observed with 
IRS (e.g. NGC 4552; Smith et al. 2007). Negligible silicate 9.7/im absorption is seen in all 



cases. 



Overall the spectra are typical of normal starburst galaxies (i.e. iBrandl et al.ll2006l ). 
However one unusual feature of these spectra is the dominance of the 11.3/im PAH feature 
over the 6.2/im PAH feature in the k+a and old k+a types. To emphasis this point in Figure 
[2] we show the spectra stacked by type, and also the stacked & continuum subtracted IRS 
spectra. 

For each object key optical emission line fluxes (H/3, [OIII]5007, Ha and [Nil]) are 
determined by fitting the residua l emission after re moval of best fit stellar population models 
using the GANDALF IDL code JSarzi et al.ll2006h . 



Equivalent widths (EWs) and flu xes for the PAH fe atures in the IRS spectra are mea- 
sured via use of the IDL PAHFIT too KlSmith et al. 1120071 ). As the mid-IR coverage is limited 
to 5-12/im the determination of the continuum around 7-12/im is troublesome due to the 
degeneracy between silicate absorption and PAH emission. However as our spectra appear to 
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Fig. 1. — IRS SL spectra for the three sub-samples under investigation here. For comparison 
the starburst template of Brandl et al. (2006) is also shown. Spectra are normalised to 
arbitrary flux at 11.3/im for the sake of clarity. 
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have no signs of strong silicate absorption around 9.7/im we constrain the sihcate absorption 
be be zero by running PAHFIT with the /NO_EXTINCTION flag set. 

Table [T] lists the optical properties, while Table [2] lists the mid-IR line fluxes and equiv- 
alent widths. 

To identify correlations between these observed properties and the post-starburst evo- 
lution we try to fit a two component stellar population model to each galaxies spectrum. 
The models are similar to those used to select the sample; a young stellar population super- 
imposed on a homogeneously old single stellar population. The young component is taken 
to be a single burst with a exponentially declining star formation rate of r = 100 Myr. The 
old component is taken to be a homogeneously 11 Gyr old single stellar population. In both 
cases models are provided by the code of Maraston (2005). To investigate both the time since 
starburst and the burst fraction a total of 24 steps in burst time from 10 Myr to 5 Gyr, and 
16 burst fractions from 0.01 to 0.15 are considered. The affect of dust is included according 
to the model of Calzetti et al. (2000) by allowing E(B-V) to vary between 0-0.2 (A^, ^0-1). 
Finally three steps in metallicity are considered; Z = 0.5Zq , Z = Zq and Z = 2Zsun- In 
total 5760 models are considered. The best fit solution is found via fitting to the standard 
25 Lick indices as well as broadband fluxes from; Galex NUV & FUV, SDSS -u, g, r, i, z and 
Spitzer 3.6//m & 4.5//m. Where an object is undetected the upper limits are incorporated. 
The burst time and burst fraction are taken to be the best-fit values after marginalising over 
the other model parameters, la errors are taken to be region surrounding the maximum 
containing the 68% of the probability distribution. These quantities are quoted in Table [TJ 

From Table [1] it is clear that our k+a selection technique has been effective in isolating 
galaxies which have undergone recent starbursts. The time since burst of the young k-|-a 
sample is ~ 0.2 Gyr, the k-|-a sample ~ 0.6 Gyr and the old k+a sample ~ 1-3 Gyr. These 
values compare well to our predicted time scales of 0.03, 0.3 and 1 Gyr for the 3 classes. Also 
shown in Table [1] is the flux measured in the key AGN diagnostic lines H/?, [OIII] ,Ha and 
[Nil]. By applying AGN diagnostic criteria presented in Stasinska et al. (2006), a variation 
on the earlier work of Baldwin, Phillips & Terlevich (1981) , to these we can determine which 
process dominates the optical emission lines. A clear distinction can be seen between the 
young k+a galaxies, which are dominated by star formation, and the older classes which are 
either AGN dominated, or mixtures of AGN and low-level star formation. 
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Fig. 2. — Stacked (top) and continuum subtracted (bottom) IRS SL spectra for the three 
sub-samples under investigation here. Oka 1 is omitted from the old k+a stack as it is 
clearly AGN dominated. For comparison the starburst template of Brandl et al. (2006) is 
also shown. Spectra are normalised to arbitrary flux in the 11.3/xm line for the sake of clarity. 



Table 1. Measured optical properties of sample 



ID 


Name 


z 


Burst Time (Gyr) 


Burst strength 




[OIII]5007 


Ha 


[Nil] 


Type'' 


Young k-l-a 


Ykal 
Yka2 
Yka3 
Yka4 


J103254. 44+581638. 5 
J104826. 85+564026. 8 
J105326. 83+580956. 5 
J163332. 66+404827.0 


0.08 
0.07 
0.13 
0.08 


"■^'0.07 

"■^■^-0.03 
„ 27+0.08 
"■^'-0.07 

16 + ""* 


„ + U.33 
"■^-0.28 

5+"-l 
"■5-0.3 

"■^^-0.16 
25+°-^^ 

"■^5_(,_jg 


106.5 
141.95 
229.16 
187.53 


104.6 
146.66 
156.57 
350.97 


291.84 
450.88 
770.26 
657.17 


54.29 
92.33 
234.36 
102.81 


SF 
SF 
SF 
SF 


k+a 


Ka2 
KaS 

Ka4 


J104953. 31+564551. 3 
J110214. 89+574630. 6 
J104319. 20+563151.0 


0.182 

0.14 

0.21 


74+"-^^ 

35+°-l 

58+"-^'* 
"■5»-0.13 


15+""'* 
"■^°-0.03 

"■"=-0.04 
„ „„ + 0.26 

"■OS-o.os 


2.56 

358.86 

6.83 


23.98 

102.05 

23.99 


17.16 

1015.24 

197.1 


44.74 

430.44 

135.08 


AGN 

HYBRID 

AGN 


Old k+a 


Okal 
Oka2 
Oka3 
Oka4 


J163252. 97+402056.1 
J104314. 82+570941. 2 
J104650. 43+575923. 9 
J104834. 84+583311. 3 


0.07 
0.12 
0.18 
0.13 


0.95±r, 

74+"-2 
"■ '^0.18 

1 9 + 0.3 


02+"-^'^ 
"■"■^-0.04 

01+°-°* 
"■"^-0.01 
1+0.07 
"■ -0.08 

"■■^0.28 


19.07 
22.75 
29.5 
158.3 


22.31 
15.98 
39.05 
49.25 


15.69 
87.8 
140.06 
519.12 


30.13 
48.22 
85.88 
249.84 


AGN 
HYBRID 
HYBRID 
HYBRID 



^Flux in units of 10 "'"'^ erg cm ^ s 
'^As defined by Stasinska et al. 2006 
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4. Discussion 

The unusual nature of the 11.3/im to 6.2/im PAH feature in the more evolved k+a 
galaxies is an interesting discovery. To quantify the nature of this Figure [3] shows the ratio 
of the 6.2/im flux to the 11.3/im and 7.7/im lines as a function of time since starburst. 

The S11.3/S6.2 and S7.7/S6.2 PAH ratios appear to evolve from the values typical for 
local starbursts (Brandl et al. 2006) to higher values with increasing time since starburst. 
Interestingly an equivalent rise is not seen in the S11.3/S7.7 ratio. 

To quantify these trends we find the best fit log-linear relationship between the variables. 
As the errors in both burst time and the flux ratios are significant simple linear regression 
techniques are not sufficient; hence best fit parameters and their error are determined via a 
Monte Carlo simulation of the dataset. In each realisation the data points are "scattered" 
by the quoted error multiplied by a normally distributed random number. The analysis here 
utilises 10,000 realisations. The mean of the slope and intercept across the realisations is 
then taken to be the best fit value, while the standard deviation is the la error on these 
properties. Performing this analysis gives the best fit slope (intercept) for the S11.3/S6.2 
ratio vs. logio(Burst Time) of 1.2 ± 0.7 (2.1 ± 0.3), 1.6 ± 1.1 (4.0 ± 0.5) for the S7.7/S6.2, 
0.14 ± 0.14 (0.5 ± 0.06) for the S11.3/S7.7, 0.2 ± 0.3 (1.0 ± 0.1) for the S8.6/S6.2 and 0.94 ± 0.58 
(2.2 ± 0.3) for the Sn.s/Sg.e- This suggests that while the observed S11.3/S7.7 is consistent 
with no evolution, the S11.3/S6.2, S7.7/S6.2 and Sn.s/Ss.e ratios show a statistically significant 
non-zero evolution with burst time. 

Smith et al. (2007) hypothesise that variations in PAH ratios could be generated via 
either weaker radiation fields, resulting in more neutral PAHs, or preferential destruction 
of the smallest grains in the presence of a very hard radiation field. They find that the 
10 objects in their sample with the highest S11.3/S7.7 PAH ratios have clear evidence for 
LINER or Seyfert AGN activity at longer wavelengths, suggesting that hard radiation fields 
associated with AGN could be preferentially destroying the smaller grains responsible for 
the 6.2/zm and 7.7/im bands, leaving an excess of the larger 11.3/im emitting grains. 

Here we also find that all of the galaxies which have high 811.3/86.2/^111 PAH ratios have 
evidence of AGN activity in their optical spectra. Interestingly we do not find that the 
811.3/87.7 PAH ratio varies with time since starburst in the same way as the 811.3/85.2 ratio. 

However this is further complicated by the observed evolution in the 811.3/83.6 ratio, 
which is what would be expected from grain destruction, but should not be possible without 
an accompanying 811.3/87.7 ratio change. Even more challenging is the odd behaviour of the 
S8.6/S6.2 ratio, which shows no evolution despite the 87.7/85.2 showing a clear change with 
time since starburst. We put these irregularities down to underestimates in the error in the 
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Hybrid ^ 
AGN □ 






0.1 0.2 



0.5 



Time since burst (Gyr) 

Fig. 3. — Flux ratios of PAH features plotted against the time since last starburst. A trend 
of increasing S11.3/S6.2, S7.7/S6.2 and Sn.s/Sg.e with burst time can be seen. Each point is 
coded by its optically determined type; circles denote SF dominated systems, squares AGN 
dominated systems and triangles hybrid combinations of the two . The grey hatched region 
represents the range of ratios found for starburst galaxies from iBrandl et al.l (120061 ). The 
dashed line shows the best fit linear relationship between the parameters. 
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S.Qfim line flux, which is heavily blended with the larger 7.7/im feature and is thus badly 
effected by the poor estimates of the continuum; a serious possibility given that we only have 
5-12/im spectroscopy for these objects. 

However AGN are only one possible explanation for this behaviour of the PAH band 



ratios. Detailed modelling of the PAH emission models by iGalliano et al.l (120081 ) has shown 
that increasing the hardness of the interstellar radiation field (ISRF) will result in a sharp 
decrease in both the S11.3/S7.7 ratio and the S7.7/S6.2 ratio while destroying the smallest 
grains (< lO'^A) will have the inverse affect. In addition increasing the ionization level of 
the PAH emitting grains will decrease the S11.3/S7.7 ratio while leaving the S7.7/S6.2 ratio 
unaffected. 

Thus the results found here can be explained via the following; 



1. The ionization level is higher than normal in these objects, decreasing the 811.3/87,7 
ratio 

2. However the ionization level is so high that the smallest grains are being destroyed, 
increasing both the 811.3/87.7 & 87.7/85.2 (and thus the 811.3/85.2 as well). 

Thus the 811.3/87.7 ratio may be balanced out by the competing processes and appear to have 
"normal" values, while the 811.3/85.2 & 87.7/85.2 ratios show abnormal increases. While any 
source which can strongly ionize and destroy the PAH emitting grains could be invoked to 
explain this phenomena, given the evidence of low-level AGN activity in the optical spectra 
for these sources we conclude that this is the most likely cause for these unusual PAH line 
ratios. However this is purely speculative and other ionization sources such as significant 
populations of evolved stars may be capable of the levels of grain destruction and ionization 
seen here. The behaviour of the 8.6/im line ratios does not fit into this model, as we should 
see similar behaviour from the 88.6/85.2 ratio as the 87.7/85.2, which we do not. 



5. Conclusion 



We have performed Spitzer IR8 spectroscopy of 10 post starburst (k+a) galaxies at 
three distinct evolutionary stages. In each case the mid-IR spectra shows the signs of star 
forming galaxy with no obvious AGN activity. Interestingly we find that relative strength 
of the 11.3/im to 6.2/im, 11.3/im to 8.6/im and 7.7/xm to 6.2/im PAH features show some 
correlation with the time since last starburst, while the 11.3/im to 7.7/im and 83.5/85.2 PAH 
ratios do not. We conclude that this behaviour is most likely a result of underlying low- level 
AGN activity in these objects. 
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Table 2. PAH line Fluxes and Equivalent Widths 



ID 


Name 


6.2/Am 




7.7/um 




8.6^m 




ll.S/am 








Flux* 


EW 


Flux 


EW 


Flux 


EW 


Flux 


EW 


Young k+a 


Ykal 


J103254. 44+581638. 5 


1.7 ± 0.8 


1.99 


5.1 ± 0.8 


5.21 


1.4 ± 0.2 


1.26 


1.8 ± 0.8 


1.02 


Yka2 


.1104826.85+564026.8 


0.4 ± 0.3 


0.36 


1.8 ± 0.3 


1.11 


0.5 ± 0.09 


0.35 


0.7 ± 0.4 


0.56 


Yka3 


J105326. 83+580956.5 


4.7 ± 1.0 


1.92 


9.2 ± 1.2 


3.54 


2.4 ± 0.1 


1.0 


3.7 ± 0.8 


1.6 


Yka4 


J163332. 66+404827.0 


2.0 ± 0.6 


2.92 


4.2 ± 0.6 


3.94 


1.0 ± 0.2 


0.81 


1.4 ± 0.7 


0.65 


k+a 


Ka2 


J1049B3. 31+564651. 3 


0.6 ±0.2 


0.46 


2.1 ± 0.3 


1.8 


0.9 ± 0.2 


0.76 


1.1 ± 0.3 


0.95 


Ka3 


J110214.S9+574630.6 


3.6 ± 0.4 


1.2 


11.0 ± 1.4 


3.6 


2.6 ± 0.2 


0.91 


4.1 ± 0.8 


1.3 


Ka4 


.1104319.20+563151.0 


1.8 ± 0.4 


0.76 


6.0 ± 0.7 


2.5 


1.3 ± 0.2 


0.67 


3.1 ± 0.8 


1.4 


old k4 


-a 


















Okal 


.1163252.97+402056.1 


















Oka2 


.1104314.82+570941.2 


0.8 ± 0.3 


0.28 


3.6 ± 0.3 


1.7 


1.0 ± 0.08 


0.59 


2.7 ± 0.5 


1.8 


Oka3 


J104650. 43+575923. 9 


1.1 ± 0.2 


0.24 


5.2 ± 0.8 


1.5 


1.0 ± 0.2 


0.32 


2.2 ± 0.6 


0.89 


Oka4 


J104834.84+583311.3 


1.1 ± 0.3 


0.30 


6.3 ± 0.9 


2.1 


0.9 ± 0.2 


0.37 


2.8 ± 0.5 


1.2 



